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ABSTRACT 

Mass modelling of early-type systems is a thorny issue; even for the few close by 
galaxies for which kinematic data is available, the implementation of this data can 
get embroiled in problems that are hard to overcome, unless more complete data 
sets are available. The mass-anisotropy degeneracy is a typical example of this. In 
this paper, we present a new mass modelling formalism for ellipticals that invokes 
no other observations other than the central velocity dispersion (cto) and photometry. 
The essence of the method lies in choosing a local mass-to-light ratio {M / L) profile 
for a galaxy, with which the deprojected luminosity density distribution (along the 
major axis coordinate x) is scaled. The resulting discontinuous mass density profile 
is then smoothed, according to a laid out prescription; the local M/L profile that 
stems from this smoothed mass density, is found to be significantly different from the 
raw M/L distribution. A suite of model galaxies (both Sersic and cored in nature) 
is used for intensive experimentation in order to characterise this raw M/L profile 
and in each case, the mass density recovered from this mass modelling technique is 
compared to the known mass distribution. We opt to work with a raw M/L profile 
that is a simple two-stepped function of x, with a low inner and higher outer value of 
M/L — Tin and To„i, respectively. The only constraint that we have on this profile 
is in the centre of the galaxy, via ctq- This value of co is implemented in the virial 
theorem to obtain an estimate of the central M/L ratio of the galaxy. The fallibility 
of the virial mass estimate is taken care of, by allowing for a range in the values of 
Tin that can be used for a given galaxy model. Moreover, our experiments indicate 
that T out is uniquely known, for a given Tin] for cored galaxies, this functional form is 
found uniquely dependant on the core radius. The physical basis for such a connection 
to exist between the inner and out M/ L is discussed. The jump radius of the raw M / L 
profile is chosen to be about thrice the effective radius of the galaxy. In this way, the 
local M/L distribution is completely specified and mass profiles of ellipticals can be 
constructed till about 3 times the effective radius. The proposed technique is extended 
to predict a mass profile for M87 which is then successfully compared to distributions 
reported earlier from kinematic analyses. 

Key words: galaxies: kinematics and dynamics - galaxies: structure - celestial me- 
chanics: stellar dynamics 



1 INTRODUCTION 

The recovery of the total mass distribution in early type 
systems is a knotty problem. The obvious reason for this is 
lack of information; photometry offers access to the stellar 
mass but the deduction of the arrangement of the total mass 
in a system remains elusive, unless kinematic information is 
invoked. Even when velocity information is available, the 
implementation of the same is besotted with problems of 
varying degrees of subtleness (Merritt 1993, Merritt 1996, 
Chakrabarty & Saha 2001, Genzel et al. 2000). 
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Of the many debacles that plague the processing of 
kinematics, the worst is perhaps the mass anisotropy de- 
generacy; even in the simplistic case of a spherical system, 
the mass distribution can be known only as a function of 
the anisotropy in velocity space (Binney & Tremaine 1987). 
Getting an independent handle on anisotropy is difficult, if 
not impossible in most cases. Cote et al. (2001) attempt this 
for the system of globular clusters in M87, by constraining 
the total mass from the work of McLaughlin (1999). Their 
conclusion indicates that the velocity distribution of this sys- 
tem in M87 is "almost perfectly isotropic" . This is consistent 
with the appraisal of isotropy in globular cluster systems in 
other galaxies (Zepf et al. 2000, Cote et al. 2003). 
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An attempt is also underway to produce a fully 
anisotropic algorithm, that is similar in flavour to the exist- 
ing code CHASSIS that works by implementing kinematic 
data of tracers to recover the potential and distribution func- 
tion from which the tracer sample is most likely to have been 
drawn (Chakrabarty & Saha 2001, Chakrabarty & Porte- 
gies Zwart 2004). However, a large number of kinematic ob- 
servations is required by this proposed extension. Nonethe- 
less, what is heartening is that such ambitious data sets are 
less of a figment of the astronomer's reverie today, owing to 
the huge kinematic surveys that are planned or already in 
operation. 

Another technique that is often employed in mass mod- 
elling of elliptical galaxies involves the utilisation of X-ray 
measurements. The procedure suffers from uncertainties in 
the observed temperature profile and simplification in terms 
of the assumption of a single component in the X-ray emit- 
ting gas and of hydrostatic equilibrium. A comprehensive 
review of various the shortcomings of the conventional mass 
modelling techniques is presented in Mamon & Lokas (2005). 

Thus, it is clear that at the present moment, utilisa- 
tion of kinematic information in mass modelling, suffers from 
egregious drawbacks. Even more crucially, such data is hard 
to come by in many external galaxies, including those at 
higher redshifts. For this latter class, the implementation of 
velocity information in formalisms that stem from a stellar 
distribution function, is all the more debatable since systems 
outside the local universe are most probably not relaxed. 

This indicates an opening for a technique that is able 
to extricate the mass profile of early type systems from 
photometry, using as little velocity information as possi- 
ble. Moreover, the used kinematic information should not 
be used in a procedure that relics on the rigorous employ- 
ment of the equilibrium stellar distribution function. Such a 
scheme is advanced in this short report. 

It merits mention that any such mass modelling formal- 
ism for elliptical galaxies will bear the potential of shedding 
light on the connection between the luminous and dark mat- 
ter content in such systems. In the past, such correlations 
have been explored in spiral galaxies (Sancisi 2004, Verheijen 
2001 and Noordcrmcr 2006), but very little is known about 
the distribution of the total mass in ellipticals; Padmanab- 
han et al. (2004) and Lintott, Ferreras & Lahav (2006) dis- 
cuss sample galaxies in this context. The scheme discussed 
here will indicate the same in model systems. 

The paper is arranged as follows. The next section pro- 
vides an exposure to the details of the proposed scheme. In 
Section 3, the Sersic model considered in the paper is intro- 
duced. Section 4 is devoted to the results obtained with the 
toy galaxy models, while the following section (Section 5) 
delineates the galaxy models that can be tackled with the 
proposed scheme from those that are outside the purview 
of this scheme. Section 6 discusses the details of the results 
obtained with the data for M87, as obtained from the ACS 
Virgo Cluster Survey (ACSVCS), (Ferrarese et al. 2006), 
and comparison of the recovered mass model with those pre- 
dicted by Cote et al. (2001) and Romanowsky & Kochanek 
(2001), (from the analysis of kinematic data of tracer popu- 
lations in M87) . The last section is devoted to the discussions 
of some of the salient points that were raised in the other 
sections and a summary of the results. 



2 DETAILS OF THE FORMALISM 

We work under the assumption that the surface brightness 
profile of a galaxy is known to us. The central velocity dis- 
persion ((Jo) is also considered known from observations. 

2.1 Deprojection 

The first step is to deproject the observed brightness pro- 
file into the three dimensional luminosity density distribu- 
tion. This can be done with the non-parametric deprojec- 
tion algorithm DOPING, which is the acronym for Depro- 
jection of Observed Photometry using an INverse Gambit 
(Chakrabarty & Ferrarese 2006). 

This algorithm can perform the deprojection under a 
general geometry that needs to be specified as an input. 
The other fundamental input is the inclination. The mear- 
surements that DOPING processes are the surface bright- 
ness measurements and the projected ellipticity profile, at 
different values of the major axis coordinate x, i.e. the sur- 
face brightness map of the galaxy. This is used to provide 
the three dimensional luminosity density distribution, the 
line-of-sight projections of which arc the best reproductions 
of the brightness data, along any azimuth on the plane of 
the sky. The deprojection of the observed photometry for 
systems with varying intrinsic shape requires the implemen- 
tation of a regularisation technique, though systems in which 
the ellipticity is uniform are more easily dealt with. It needs 
to be mentioned that DOPING can in principle, deal with 
changes of eccentricity and twist angle, with galactocentric 
distance. 

The errors on the deprojected luminosity profile are the 
errors of the analysis, which measures the zbl-u extent of 
the excursion of the algorithm around the global maxima in 
the likelihood function; likelihood is maximised on spotting 
the luminosity density distribution that projects best into 
the data. The observational errors can be also incorporated 
into the analysis. 

2.2 Models 

In this paper, we attempt to mass model a set of toy elliptical 
galaxies. In terms of measurable quantities, these model sys- 
tems are ascribed a brightness profile I{x) and a projected 
ellipticity e, that we choose to fix at 0.4, (implying that for 
edge-on viewing, the intrinsic eccentricity e=0.8). All the 
toy galaxies are assumed to have been observed in the SDSS 
7-band and placed at the distance of Virgo (considered as 
17Mpc). 

Multiple runs are carried out with the toy galaxies being 
assigned varying values of central velocity dispersion (uo), 
from which the local mass-to-light ratio at the centre of the 
galaxy is calculated in a way discussed below. These observ- 
ables are the inputs into the mass modelling system that 
is proposed in the paper. The recovered total mass density 
distribution is then compared to the true mass density pro- 
files {pt{x)) of these models. The degree of overlap between 
the predicted and known mass distributions is then visually 
tracked, in order to confirm the success and limitations of 
the proposed mass modelling scheme. 

The models are described by a luminous component (lu- 
minosity density L{x)) lying embedded in a dark halo that is 



Mass Modelling with Minimum Kinematic Information 3 



chosen to be NFW type (density: pnfw{x). Thus, the true 
total mass density of a model galaxy can be written as: 



1000 P 



pt{x) = pNFw{x) + alQ 2.5 ^ 
where the dark matter density is 

pNFW{X) — 



(1) 



(2) 



and a is a measure of the luminous matter fraction, such that 
its reciprocal indicates the halo fraction in the toy galaxy. 
Also, 7/0 is the absolute magnitude of the Sun in the same 
waveband as the brightness observations. 

We treat two sets of model brightness distributions: the 
cored and the Sersic type profiles. As far as the cored bright- 
ness distributions are concerned, the following prescription is 
used. The brightness profile of the test galaxies is extracted 
from the line-of-sight integration of an analytically chosen 
luminosity density function. (Thus, the deprojection of this 
brightness data can be compared to the known luminosity 
density to confirm the correctness of the deprojection). This 
analytical luminosity density distribution Lt (x) is chosen to 
have the form: 



Lt{x) 



A 



[r-2 + 2:2 + yV(l-e2)+z2]' 



(3) 



where A is the central luminosity density and the core ra- 
dius of the test galaxy in question. Thus, various toy surface 
brightness data sets are prepared, with varying amplitude 
and core radii. The deprojection is carried out under the 
assumptions of oblateness and an edge-on viewing angle. 

We normalise the deprojected luminosity density pro- 
files of all the model galaxies such that the central value of 
L{x) is about 100 Lqpc"'^, leaving the models distinguish- 
able from each other in terms of the core radius only. Each 
such model luminosity density profile is then ported to the 
mass modelling formalism (discussed below) and the out- 
come is compared to the known mass density distribution 
of this toy galaxy. Figure 1 shows the normalised luminosity 
distributions of the cored galaxies (of core radius 7.07pc and 
about 22.4pc) that have been used in our experiments. 

For the Sersic models, the brightness is chosen to be de- 
scribed by a Sersic index of 5.3 and a scale length of about 
108pc. This central density is normalised to 1000 L0pc~^, 
i.e. 1 order of magnitude higher than the normalisation fac- 
tor for the cored galaxies. 



2.3 Raw Two-stepped Discontinuous M/L Profile 

Once the deprojection is performed, the next stage is to 
choose a total local mass-to-light ratio (M/L) profile that 
scales the obtained luminosity density profile in such a way 
that it would hopefully emulate the total mass density. Or 
at least, such is the aim, as in other procedures that hope to 
connect photometry to a mass model. We suggest selection 
of a discontinuous, two-step local M/L: 



M/L = Tin ifa; < Xi„ 
M/L = Tout ifa; > Xin 



(4) 



The luminosity density is scaled by this mass-to-light ratio 
profile and then smoothed. Of course, the success of the 
scheme depends on the selection of Xin, Tin and Tout. 




100 1000 
DisLance along x (pc) 

Figure 1. Deprojected luminosity density profiles of toy galaxies 
with core radii of about 7.1pc (in black) and 22.4pc (in magenta), 
recovered by DOPING from model brightness profiles that were 
constructed from the analytical luminosity distribution given in 
Equation 3, normalised to 100 Lqpc"^. The Sersic luminosity 
density is shown in cyan, normalised to 1000 L0pc~'^. 



Additionally, we need to be cautious that this method is 
susceptible to spurious answers at distances where the dark 
matter halo dominates the system. Thus, we aim to produce 
a mass profile over a radial range that is not too extensive 
(not exceed about 3Xe, where Xe is a measure of the effec- 
tive radius of the galaxy). We ascertain Xe by measuring 
the slope (m) of the straight line that is fit to the plot of 
core-removed logio(I) against x^^^: Xe = (— 3.33/m)''. 

So, what would amount to a judicious prescription for 
the selection of the free parameters in this formalism? If we 
were to assume that mass follows light in elliptical galaxies, 
then we could set Tt„ — Tout and Xin could be set to any 
value in the range of [0, 3Xe]. But there is no reason to 
believe that mass does follow light in elliptical galaxies. In 
fact, gas kinematics and X-ray studies indicate otherwise; 
Buote et al. (2002) report that the hypothesis that mass 
follows light can be rejected at 96% for their oblate models 
for NGC720 and 98% for their prolate models. Bertola et al. 
(1993) point out that the run of the mass-to-light ratio with 
galactocentric radius in ellipticals is similar to that in spirals. 

Kinematics too rule out the possibility of a constant 
M/L at 95% confidence (Gerhard et al. 2001). More re- 
cently, the analysis of the kinematics of planetary nebulae 
observed in five elliptical galaxies, by the Planetary Neb- 
ula Spectrograph, has led to the announcement that these 
galaxies are drastically short of dark matter in the outer 
parts (Romanowsky et al. 2003, Romanowsky et al. 2004). 
Napolitano et al. (2005) advance a parameter that measures 
the logarithmic slope of the mass-to-light ratio in a sample 
of early type galaxies, to find that it increases with the stel- 
lar mass. They conclude that the distribution of dark matter 
or its fraction changes with radius. 
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2.3.1 Jump Radius 

Thus, we are left to make as good a choice for the parameters 
as possible, from whatever can be learnt about the system 
from the minimum of observables. Now, there is no obser- 
vational constraint that can suggest a value for the jump 
radius Xin, at which the above two-step local M/L profile 
should suffer its discontinuity. However, one obvious length 
scale in the problem is Xe. So we can choose to impose a 
certain multiple of Xe as the jump radius; in fact, experi- 
ments indicate that a good value of Xin is 3Xe. This is not 
a particularly bad choice cither, given that by 3 effective 
radii, the dark matter content can be expected to begin to 
exert its influence. Having fixed Xin, we now proceed to use 
whatever observable information there may be available to 
constrain Tj„ and Tout. 

2.3.2 Inner M/L 

One measure that is often available from observations is the 
central velocity dispersion in galaxies, ao can be utilised to 
get a handle on the local mass-to-light ratio in the inner 
parts of the galaxy (i,e, Ti„). This can be done by im- 
plementing (To that is measured aX, x = xo (say), in the 
virial theorem. However, the estimation of the mass enclosed 
within xo {Mq), from which the average mass density (po) 
inside xq can be calculated in this way, will be mired in 
errors owing to a lack of information about the degree of 
anisotropy present at x = xo, as well as the extent of de- 
viation from sphericity at xq. Virial theorem suggests that 
under the assumption of velocity isotropy, 

„ 2 GMo 

ScTo = — (5) 

and oblateness and edge-on inclination imply 

Mo = ^po4 Vl - (6) 

This estimated po will deviate from the true inner mass den- 
sity, as explained above. Anisotropy will spuriously lower the 
estimated enclosed mass and therefore the calculated po, for 

a given ao. For example, Padmanabhan et al. (2004) sug- 
gest that the dynamical mass enclosed within the projected 
half-light radius {Re) of the massive ellipticals in their SDSS 
sample is given as 1.65^ times aQR^/G, with a systematic 
error of 30% (Lintott, Ferreras & Lahav 2006). 

Expectedly, the central M/L {M/L{0)) has a one-to-one 
and monotonically incrcEising relationship with the luminous 
matter fraction index a (lower left panel of Figure 2). As 
shown in the figure, M/L(Q) is well represented by q, though 
the deviation of [M / L{Q)]/ a from unity exists - this ratio 
decreases from about 0.22 to 0.07 as we increase from 
7.1pc to 22.4pc. We will use this result to justify our choice 
of a for parametrising the estimated central M/L. 

We scan through a range of estimated central M/L val- 
ues (i.e. a) and monitor the success of the proposed tech- 
nique; in the process, the ranges of Tin and Tout correspond- 
ing to a given choice of the central M/L value are recorded. 

Now, experiments with our model galax;ies of a given 
core radius Vc indicate that for a given a, there exists a 
range of values for Ti„ that will result in compatibility of the 
estimated and known mass density profiles (upper left panel 
in Figure 2). This band of allowed Ti„ values for a given a, 



varies with core radius of the normalised luminosity density 
profile in hand. For our cored galaxies, this band of Ti„ 
values that are allowed at a given a, range from just above 
the corresponding value of a to approximately 2 times this 
a, with this factor varying only weakly with core radius Vc. 

If we assume the value of 2.72 between the true value 
of the mass enclosed within an effective radius and a^Re/G, 
as suggested by Padmanabhan et al. (2004) (and used by 
Lintott, Ferreras & Lahav 2006), then the true central M/L 
should be about 1.82 times the virial estimate. Thus, the 
range of 2a that is allowed in Ti„, for a given estimate of 
M/L{0), appears to be inclusive of real ellipticals, at least 
the systems that are as massive as those in the sample of 
Padmanabhan et al. (2004). 

If it so happens that a certain galaxy is far too 
anisotropic at the centre, such that the true M/L{0) ex- 
ceeds the upper bound on the Ti„ allowed at the a at hand, 
then the mass density profile that will result from using the 
highest value of the permitted Ti„ in our mass modelling 
formalism, will not concur with the profile that results from 
using the lowest allowed Ti„. Thus, in such circumstances, 
when our scheme fails, we will be aware that it has done so. 

Thus, it may be summarised that oven though the virial 
theorem is invoked to help constrain Ti„, there is enough 
scope available within the formalism to accommodate for 
the problems caused by velocity anisotropy in traditional 
kinematic modelling. This and other flexibilities offered by 
the proposed formalism are discussed in greater detail in 
Section 7. 

2.3.3 Outer M/L 

Our experiments indicate that it is possible to achieve a 
one-to-one correlation between Ti„ and Tout, for a toy con- 
figuration of a given core radius. In other words, irrespective 
of a, whenever a value of Tin is chosen. Tout is fixed. This 
relationship between Tin and Ti„ is depicted on the top 
right panel in Figure 2. 

Actually, there exists a range of Tout for a given Ti„, 
with the extent of this range maximising itself at values of 
Tin that lie in the middle of the range that is allowed for a 
given a. However, to keep matters simple, we present only 
the value of Tout that lies at an extrema of this range, for a 
given Tin. 

Analytical fits to the plots suggest that Tout varies as 
Al exp[— Tj„/r] + Ao, where Ai and arc constants that 
are determined from fitting this exponentially decaying func- 
tion to the plot of Tout against Tf„; r measures the decay 
rate. Fitting this functional form to the experimentally ob- 
tained plots, suggest power-law dependences of both the ad- 
ditive and multiplicative constants, on re. 

(i) Ai^re-'-'\ 
(u) Ao ~ rc"^ '^® and 

(iii) T has a weak dependence on rc, according to r = 
l.OOS-'^r;* + 1.85. 

That the inner and outer local M/L ratios are con- 
nected, betrays a relationship between the luminous and 
dark matter content in these toy elliptical galaxies. This 
facet of the models and the conclusions that we can draw 
on the basis of these results for elliptical galaxies in general, 
are further explored in Section 7. 
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Figure 2. Lower left panel: the luminosity matter fraction a, plotted as a function of the estimated central M/L. Results for rc=7.1pc 
are shown in black, for 14.4pc in red, for 17.3pc in green, for 20pc in blue and for 22.4pc in magenta. Top left panel; figure showing the 
allowed ranges in Ti„, as a function of a (which measures the estimated central local M/L), for cored galaxies of varying core radii; the 
colour indexing used here is the same as above. Top right panel: Tout plotted as a function of Ti„. 



2.3.4 Smoothing 

Once the local A'l/L is specified, it is used to scale the lumi- 
nosity density distribution, in order to obtain a mass density 
distribution, which is indeed rendered discontinuous by the 
process of its formulation. Thus, the next stage is to smooth 
this mass density distribution. 

We choose to smooth this density profile with a tri- 
angle filter of filter size given by 2Xe, i.e. two successive 
applications of a box filter of filter size corresponding to the 
length Xs. It is appropriate that the smoothing filter size 
be a function of one of the characteristic scale- lengths of the 
system; Xe is the most obvious candidate that satisfies this 
criterion. Thus, the trends in Tout with variation in Ti„, 
for the different normalised toy galaxy configurations, have 



been extracted from experiments performed with a smooth- 
ing filter size of Xe. However, other smoothing windows are 
also usable - except that the exact nature of the relationship 
between Tout and Tout will then be different from what is 
prescribed above, for a given Tc. 



3 SERSIC GALAXIES 

All the toy models that are discussed above correspond to 
examples of surface brightness profiles that are cored. In or- 
der to investigate the full applicability of our mass modelling 
technique, a Sersic type galaxy is now examined. 

We consider a toy Sersic surface brightness distribu- 
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tion, described by a central density of 1000 Lqpc~^ , a scale 
length of about lOOpc and a Sersic index of 5.3. When de- 
projected under the assumption of sphericity, a luminosity 
density profile along the semi-major axis is recovered. The 
central density chosen for this toy model is 1000 Lqpc"'^ 
instead of 100 Lqpc"^ (as in the cored models), since the 
proposed mass modelling technique does not work for Sersic 
models at the lower normalisation. 

Once normalised in this way, a range of central (local) 
mass-to-light ratio values are imposed, in order to scan the 
range of parameters over which our scheme is viable. As 
before, the central M/L values are imposed in terms of the 
parameter a (the luminous matter fraction). The raw two- 
step local Af/L profile that is imposed upon this luminosity 
density is again chosen to be described by an inner M/L of 
Tin and an outer M/L of Tout, with the jump radius given 
by 3Xe. 

In these Sersic models, Xe is defined as the slope (m) 
of the straight line that is fit to the plot of logio(I) against 
x^^^: Xe = (— 3.33/m)'*. For the toy model at hand, we find 
Xe ?»726.8pc. 

Once the two-stepped local AI/L profile is fully char- 
acterised, the deprojected luminosity density is scaled by it 
and the resulting mass density distribution is smoothed. As 
before, it is smoothed by two successive applications of a 
boxcar smoothing filter of filter size given by Xe. 

For a given choice of the central M/L, a range of Tin 
are allowed. This is shown on the left panel of Figure 3. 
This range is evidently smaller for the Sersic toy model, 
compared to the cored models. Again, we find that a one- 
to-one relationship exists between Tin and Tout, as shown 
in the right panel of Figure 3. 

Thus, it is apparent that galaxies with brightness distri- 
butions that are best described by a Sersic profile, can also 
be tackled by the proposed mass modelling scheme. The two 
main differences in the formalism, as applicable to the cored 
and Sersic brightness profiles are the following. 

(i) The luminosity density for a Sersic-type galaxy needs 
to be normalised to 1000 Lqpc"^, while a cored brightness 
profile would require a normalisation factor that is 1 order 
of magnitude smaller. 

(ii) For the determination of X^ of the Sersic-type galax- 
ies, a straight line fit to the plot of logio(I) against x^^'^ 
needs to be performed while for the cored galaxies, a linear 
fit is sought to the relationship between the core-removed 
logio(I) profile and x^^"*. 



4 RESULTS 

Figure 4 brings out the change in the local mass-to-light 
ratio profile, as a result of the smoothing; the mass density 
distribution that results from the smoothing of the scaled lu- 
minosity density profile, is compared to the luminosity den- 
sity, to extract the smoothed local M/L. This is seen in the 
figure to be different from the raw, two-step, discontinuous 
M/L distribution that is used to scale L{x). The smoothing 
is noted to afi^ect mass-to-light on both sides of the jump 
radius. 

Figure 5 displays the degree of overlap between the esti- 
mated and true mass density profiles, for toy cored galaxies 
that are normalised to the same central luminosity density 




1000 



Figure 4. Plot of the raw (in broken red lines) two-stepped lo- 
cal A// L ratio that is used to scale the luminosity density profile 
of the toy galaxy of core radius 17.3pc, given model parameters 
of a=3, Ti„=6 and Tout=31. The resulting discontinuous mass 
density profile is then smoothed according smoothing prescription 
mentioned above. When this smoothed mass density distribution 
is compared to the deprojected luminosity density profile, the re- 
sulting local M/L profile is shown in filled circles. As is apparent, 
the raw two-stepped M/L distribution varies significantly (be- 
yond error bars) from the final smoothed M/L ratio distribution. 



and characterised by varying shapes, i.e. varying rc (upper 
panels). The lower panels of the same figure indicate a simi- 
lar comparison, when a is made to change. In all these cases, 
when we utilise the prescribed values of Tin and Tout, the 
predicted mass density profile reproduces the true mass den- 
sity profile very well, within error bars. 

Figure 6 brings out the resemblance between the known 
mass density distribution along the major axis (in green) 
and its predicted counterpart (in red), for the Sersic model 
investigated, for a value of Tin that is picked from the range 
allowed for a given choice of the estimated central M/L. 



5 ALLOWED GALAXIES 

This scheme of characterising the two-step local mass-to- 
light ratio profile is found to work only as long as the galaxies 
at hand lie within a certain band in the Ms—rs space, where 
Ms is the mass scale and is the scale radius that describe 
the NFW dark halo mass density profile (see Equation 2). 
This band is translated into a range of the allowed values 
for the cored galaxies, in the V200 — rs space, (where V200 is 
the circular velocity at radius r2oo), in Figure 7. The broken 
line in this figure indicates the state of the NFW halos, as 
given in Hoekstra, Yee & Gladders (2004) . The range in Ms 
allowed for Sersic-type galaxies is slightly more constricted, 
as compared to the cored galaxies. 

To clarify, by these "allowed ranges", we imply the 
ranges in V200 and rs, over which the same Tin and Tout 
values are valid. 
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Figure 3. Left panel: the allowed ranges in Ti„, as a function of the estimated central local M/L, as parametrised by a, for the Sersic 
model that we consider. Top right panel: Tout plotted as a function of Tj„; the solid line represents an exponentially decaying function 
of Tin that is a very good fit to our experimentally obtained data. 



Thus, it is apparent that the mass modeUing scheme 
that is being proposed here is applicable even to systems 
that may exhibit a significant range in concentration. At 
the same time, we note that dwarf galaxies which typically 
have ^200 ~ 100 kms~^ cannot be characterised by this for- 
malism. 

In addition to these constraints on the Ms and Vs, the 
gala:xies also need to be "ordinary" in the sense that they 
should, qualitatively speaking, have a lower inner mass-to- 
light ratio than a higher M/L at higher radii. However, if 
the galaxy at hand is known to have a central mass con- 
densation from independent measurements (such as M87), 
the raw two-step M/L profiles discussed above (described 
by a low Ti„ and a Tout) will not be apt. We refer to such 
galaxies as "not ordinary" in this paper. For such a system, 
an alternative three-step, local mass-to-light profile can be 
conjured. In such a configuration, the value of the M/L in 
the innermost region must be retained as high, while the 
intermediate range is envisaged to be described by a low 
M/L. Outside of the second jump radius, the dark matter 
contribution is expected to kick in; therefore, the mass-to- 
light ratio in the outermost part should be held at a high 
value. The values of the M/L in the different parts of such 
systems are not expected to be given according to the pre- 
scriptions mentioned above; nor arc the two relevant jump 
radii in these cases, expected to be given as multiples of Xg 
in the same way as for "ordinary" galaxies. 

In general, the mass modelling of such systems axe be- 
yond the scope of the proposed formalism. Nevertheless, 
in Section 6, we will show how the basic concept behind 
our proposed mass modelling formalism can be extended to 
model M87, using whatever we know of the physics of M87, 
from literature. 



6 M87 

In this section, we present a mass model of the galaxy M87. 
Wc compare our results with the enclosed mass profile put 
forward by Cote et al. (2001) and the mass density of M87, 
suggested by Romanowsky & Kochanek (2001). 

M87 (NGC 4486) is a cD galaxy that lies near the cen- 
tre of the Virgo cluster. It heis been extensively studied 
both for the photometry of its stars and its globular cluster 
system (Strom et al. 1981) as well as for their kinematics; 
some of the latest investigations are due to Romanowsky & 
Kochanek (2001), Cote et al. (2001). M87 is marked by the 
existence of a central supermassive black hole (Macchetto 
et al. 1997, van der Marcl 1994) and a dearth of dark mat- 
ter in the outer parts; we configure the details of the raw 
M/L profile accordingly. 

Given that M87 is "not ordinary" , according to the pre- 
scription suggested in the previous section, it is expected to 
require a three step discontinuous raw M/ L profile, in which 
the local M/L in the innermost step is high (in conjunction 
with the reported central supermassive black hole) and that 
in the intermediate and outermost steps are expected to be 
moderate and high respectively. However, given that M87 
has been modelled with no dark matter halo in the past 
(McLaughlin 1999), with most of its dark matter being that 
of Virgo, it makes sense for us to maintain the same value 
of the local M/L in the outer two steps (at least upto the 
radius to which photometric information is available). Thus, 
the raw M/L profile is rendered two stepped and we are re- 
quired to judiciously choose the jump radius Xin that marks 
the end of the regime where the influence of the central black 
hole dominates. The details of the choice of this raw M/L 
profile are discussed below: 

* IS chosen to be 0.2 ; the ACS surface brightness pro- 
file of M87 indicates a steepening of the profile from about 
0.2 inwards. 
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Figure 5. Comparison of recovered (in red) and model (in green) mass density profiles for the cored toy galaxies with varying a (lower 
panels) and core radii (upper panels). The two different core radii corresponding to the plots on the upper right and left are about 
17.3pc and 7pc, respectively. These toy configurations have been ascribed the same a of 10. The lower right and left panels correspond 
to the same rc of 20pc and a=3 and 7, respectively. The range of x over which the profiles are shown corresponds to thrice the Xe for 
the rc=20pc model. This is in excess of 3Xe for the models with the lower core radii. That the overlap between the recovered and true 
profiles is notable over this whole x range, even in the cases of rc=17.3pc and 7pc, refiects the realisation that the mass modelling is 
often successful even beyond x = SX^. 



n 

• Outside the central 0.2 , the mass to hght ratio is con- 
sidered to be constant, with a moderately large value of 
about 10. This is maintained till the end of the radial range 
to which the photometric data extends (about 100 ), which 
is just slightly higher than the Xe of about 96 . We find 
that the allowed range of this intermediate M/ L is 6 to 13. 

• The observed central line-of-sight velocity dispersion of 
M87 (van der Marel 1994) is about 400kms"^ This trans- 
lates to an enclosed mass of about 1.6x10^Mq. Comparison 
with photometry implies a local central M /L of about 200. 

Now, let us discuss the predicted mass models for M87, 



against which the recovered profile is checked. Romanowsky 
& Kochanek (2001) adopt an orbit modelling approach to in- 
fer the mass density, from measurements of velocities of stars 
(to about 1.5 effective radii) and globular clusters (to about 
5 effective radii). The available surface brightness and kine- 
matic observations of the stars and globular clusters over as- 
sorted radial ranges were used in this work. Romanowsky & 
Kochanek (2001) use the Schwarzschild orbit library method 
to work with generalised mass models that include a contri- 
bution from the luminous matter; this is represented by the 
luminosity density i^*(r) that is scaled by the M/L. The 
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Figure 6. Similar to Figure 5, except that in this case, the comparison is shown for two distinct Sersic models that have been ascribed 
o of 2 (left) and 9 (right). 



dark matter density is considered to be NFW type. The fit- 
ting is conducted on the joint constraints from the stellar 
and globular cluster velocity data, in a Jeans equation for- 
malism. In this work, the availability of only one component 
of the velocity vector (the line-of-sight component), must 
render a wide range of solutions possible. 

The best fit model from Romanowsky & Kochanek 
(2001) is presented in the right panel in Figure 8 in green. 
Even this best fit model is reported to suffer from the ail- 
ment that the values of the normalisation of the M87/ Virgo 
mass profiles, as obtained from the stellar and globular data 
sets separately, are unequal. This could be explained by ei- 
ther a larger, less concentrated halo or by the relaxation of 
the assumption of sphericity. In fact, the ACSVCS image of 
M87 shows clearly that the eccentricity in M87 is not zero; 
the eccentricity is seen to dip sharply in the inner regions 
and then appears to tend towards a maximum of about 0.4, 
as radius increases. 

The dynamical analysis of the radial velocities of 278 
globular clusters of M87 is presented in Cote et al. (2001). 
The total (projected) radial range spanned by the kinematic 
data set of Cote et al. (2001) extends from about 36 to 
about 6.5 effective radii. This kinematic data is then fed 
into the Jeans equation under the assumption of a spherical 
geometry. In order to evade the problem of mass-anisotropy 
degeneracy. Cote et al. (2001) take the distribution of the 
total mass of M87 and the Virgo cluster that was deduced 
by McLaughlin (1999). This enclosed mass distribution is 
shown in black in the left panel of Figure 8. In this model, 
the dark matter is assigned wholly to the cluster. 

The surface brightness and eccentricity profile of M87, 
observed as part of the ACS Virgo Cluster Survey were used 
as inputs to DOPING which recovered a luminosity profile 
for the galaxy. When the inferred light distribution is scaled 
by the two step M/L profile discussed above, and the conse- 
quent discontinuous mass density distribution is smoothed 



200 - 



180 - 



140 - 



120 - 




1-1 16 18 ao 24 

(kpc) 



Figure 7. The shaded region indicates the section of the V200 — ''s 
plane that corresponds to galaxy configurations that can be suc- 
cessfully modelled with the proposed mass modelling scheme. The 
solid black line is the NFW prediction, according to Hoekstra, Yee 
& Gladders, (2004). The red stars mark the extremes of the ranges 
in V200 Eind Vs that are allowed for the Sersic model considered 
in the paper. 



(according to the smoothing prescription suggested in Sec- 
tion 2.3.4; smoothing filter size corresponds to Xin). The 
smoothed mass density profile is used to extract the en- 
closed mass profile, which is compared to the enclosed mass 
distribution of M87 that has been predicted by Cote et al. 
(2001) (see Figure 8). The comparison appears favourable 
indeed though the comparison of the recovered mass den- 
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Figure 8. Figure showing the enclosed mass profile (left), and mass density (right) of M87. The ACS photometry for M87 is deprojected 
by DOPING to result in the luminosity density distribution of this galaxy. The recovered luminosity density is then scaled by a raw 
two-step local M/L profile that has a jump radius of 0.2 , Ti„=200 and Tout=6. The resulting discontinuous mass density distribution 
is smoothed (two successive applications of a boxcar smoothing filter of filter size 26 which corresponds to the jump radius of 0.2 ). This 
smoothed mass density distribution (in red) is compared (in the right panel) to that predicted by Romanowsky & Kochanek (2001) (in 
green). The enclosed mass profile that is calculated from this recovered mass density profile is shown in red in the left and is compared 
to the predictions of Cote et. al (2001) (in black). 



sity to that predicted by Romanowsky & Kochanek (2001) 
is less encouraging in the central parts of the galaxy. This 
discrepancy however does not reflect badly on our mass mod- 
elling technique since the density predicted by Romanowsky 
& Kochanek (2001) is also in deficit of the enclosed mass pre- 
dictions of Cote et al. (2001), due to the lack of inclusion of 
a central black hole in their model. 



7 SUMMARY AND CONCLUSIONS 

We propose a scheme to retrieve the total mass density to 
about 3 times the effective radius in ellipticals, using noth- 
ing but photometry and the central velocity dispersion: the 
exact two measurements that are most commonly available 
for most galaxies. The basic framework of this formalism 
relies upon the construction of a two-step, discontinuous, 
local M/L profile, which is chosen to suffer its discontinu- 
ity at 3Xe. Here, Xe is the major-axis coordinate equivalent 
of the effective radius, as defined in terms of the fitting of 
the observed brightness profile to the de'Vaucouleurs i?^''* 
law. Prescriptions are offered for constraining the inner and 
outer amplitudes of this sought two-step M/L profile. While 
Tin is found to lie within an allowed range of mass-to-light 
values for a given estimate of the central mass-to-light in 
the system. Tout is fixed for a given Ti„, for a galaxy of a 
given shape. The details of the prescriptions are found to 
differ slightly, depending on whether the observed surface 
brightness is noted to be cored or Sersic in nature. 

This formalism is designed to accommodate for the 
usual problems that are associated with our lack of knowl- 
edge of the state of anisotropy in the system and deviations 



from sphericity. Thus, the measured central velocity disper- 
sion is implemented in the virial theorem to pin down the 
(local) central mass-to-fight ratio, but the possibility of mis- 
calculating the same due to anisotropy always exists. How- 
ever, this mistake can be covered for, to some extent, in 
the sense that as long as the true central mass-to-light ra- 
tio (Tin) lies within a range above the estimated central 
M/L (parametrised by a), the recovered mass density pro- 
file is found to be compatible with the true distribution. 
This range is given as 2a, the central luminosity density of 
which is normalised to 100 L0pc~^. If on the other hand, 
the brightness profile is Sersic type, the normalisation fac- 
tor is 1000 L0pc~'^ and for the only Sersic model that was 
explored, the upper bound on the allowed range for Ti„ is 
relatively more constricted, viz. 1.65 times a. A one-to-one 
relationship between Tin and Tout exists in both cases. 

Such a prescription can fail if the degree of anisotropy 
in the centre of the system is so high that the virial esti- 
mate of the central mass density is rendered different from 
the true central density by a factor in excess of 2a in the 
case of cored galaxies and l.GSa, for Sersic galaxies. The 
good news is that when this happens, the formalism can in- 
dicate the same; such configurations cannot be handled by 
the proposed scheme. However, it is unlikely that the galaxy 
will be this anisotropic at the centre since a careful analysis 
of a sample of over 2000 galaxies in SDSS indicates that the 
dynamical mass enclosed with IR^ is about 1.82 times the 
virial mass estimate, as given in Equation 5 (using the result 
of Padmanabhan et al. 2004). 

An important factor that dictates the success of our 
mass modelling scheme is smoothness of dispersion profiles 
in real galaxies. As Gebhardt et al. (2000) suggested, the 
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largest variations in normalised dispersions occur in galaxies, 
either in the inner parts (R < 0.5Xe) or at radii exceeding 
2Xe. In the radial range of 0.5Xe < R < 2Xe, the variar 
tions in dispersion profiles is reported to be less than 5%. 
Gerhard et al. (2001) suggest that the M /L profile starts to 
pick up around 0.5-2Xe. This uniform dynamical character- 
istic of real ellipticals is exploited in the production of the 
formalism. 



7.1 Core Fitting Method 

At this point, it merits mention that though one part of the 
proposed mass modelling technique resembles King's core 
fitting method (King 1966, Rood et al. 1972) in some way, 
the whole of the advanced technique is very different from 
the core fitting method. The resemblance lies in the invoking 
of the virial theorem in the extraction of the mass that lies 
enclosed within a pre-fixed radius; in our work, this is the 
minimum radius (a;o) at which a velocity dispersion measure 
is available. In the implementation of the virial theorem, 
isotropy and sphericity are assumed, as in the core fitting 
method. 

But the essential difference between the two methods 
lies in the fact that while the virial estimate of the mass 
is advanced as the answer in the core-fitting method, in 
our method, this mass is used to spot the range of allowed 
Tin values that correspond to this estimate. In other words, 
given a virial estimate of the mass within xo, we attempt to 
decipher the range of values within which the true central 
(local) M/L lies. The availability of a range in this M/L 
offers some respite from the difficulties caused by the pos- 
sible existence of anisotropy, which the core-fitting method 
is not safeguarded against; in fact, it appears from work 
done with real galaxies that this respite is realistic indeed. 
Moreover, our method would forewarn us in case the system 
is so severely anisotropic that the method would fail. The 
core-fitting method has no such inbuilt safeguard. 

Most importantly, once we have chosen our Tin, (which 
then fixes Tout), smoothing is implemented. The smoothed 
mass density distribution is then found to correspond to 
a local M/L profile that is found to be different from the 
choice of Ti„, in the inner parts of the galaxy. The whole of 
the M/L distribution, from a; = to a; = SXe, is rendered 
affected by the overall smoothing. This is clearly borne by 
the Figure 4. On the other hand, the M/L in the central 
parts of the galaxy, as found from the core-fitting method, 
would of course be just what the virial estimate indicates. 

All in all, we can say that our method does invoke 
the virial theorem to translate velocity dispersion informa- 
tion into an enclosed mass value, as is done in the core- 
fitting method, but in our technique, we then go ahead to 
refine that virial estimate in multiple ways, unlike in the 
core-fitting method. The proposed mass modelling formal- 
ism is applicable to aspherical galaxies that are at least as 
anisotropic in the centre as massive ellipticals in the SDSS, 
even when these systems bear constant surface brightness 
cores, and it provides a mass distribution over a radial range 
extending to about 3 times the effective radius. The core- 
fitting method is much more constrained in the sense that 
it provides the mass profile only in the central regions of 
galaxies that do not harbour a constant surface brightness 



core and which when aspherical and/or anisotropic, cause 
the estimated mass distribution to be erroneous. 



7.2 The Case of M87 

The fundamental motivation for the raw two-stepped local 
M/L profile to be ascribed a low inner and high outer M/L 
amplitude is to replicate the anticipated trend of the dark 
matter kicking in at large radii {iXe and higher). But if a 
system is known (from independent measurements) to har- 
bour a central mass condensation then the innermost radial 
range should bo assigned a high local M/L indeed. Such is 
the case for the galaxy M87. Additionally, M87 is known to 
be a cD galaxy lying at the centre of Virgo. Thus, it can 
be modelled as being bereft of any dark matter, with all 
the dark matter ascribed to Virgo itself (McLaughlin 1999). 
Such information about the reported state of the matter 
distribution in M87 can motivate us to model this galaxy, 
using our mass modelling scheme, with a two step M/L pro- 
file in which the inner M/L amplitude is high, compared to 
the rest of the radial range covered by photometry. Once 
we configure our raw M/L profile in this way, the recov- 
ered mass distribution in M87 is found to closely follow the 
mass models predicted for M87 from earlier studies. Thus, 
the proposed technique relies upon the inclusion of as much 
measured information about the galaxy at hand, as is pos- 
sible. 



7.3 Connection between luminous & dark matter 
densities 

Once Tin is constrained. Tout is fixed; this correlation be- 
tween Tin and Tout depends on whether the system is cored 
or Sersic. A full investigation was done with the cored galax- 
ies that were normalised to a central luminosity density of 



100 Lqpc , to suggest that 
Tout = 8.32 X 10' r," + 
9.76 X 10**'r^ cxp 



(7) 



^ in 



1.005-«r*-h 1.85 



Similar analysis performed with a toy Sersic brightness pro- 
file revealed that Tout is again related to Tin via an expo- 
nential decay though a more thorough investigation needs to 
be undertaken to appreciate the dependence of this relation 
on the Sersic parameters. 

The fact that the local mass at iX^ in our toy ellipti- 
cals is related uniquely to the mass near the centre, betrays a 
relationship between the dark matter and luminous matter 
content in these systems. Given that the total mass den- 
sity in these models has a contribution from an NFW-type 
dark halo and a times a (cored or Sersic) luminosity density 
profile, (where a is a scalar), such a relationship is not al- 
together unexpected. As shown in Figure 9, it is clear that 
in the inner regions [x < Vc) the total mass density is dom- 
inated by the shape of the underlying luminosity density 
profile, while the dark halo predominates in the outer re- 
gions (a:: ~ 3Xe)- Thus, when the luminous mass is low, the 
dark mass is high and the vice versa. This trend manifests 
itself into the exponential decay of the outer M/L value as 
the inner M/L is increased. It must be kept in mind that 
this relationship between Ti„ and Tout is in the context of 
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Figure 9. Figure showing trends in dark matter density (black) 
and luminous matter density (in red) for a model cored galaxy 
with a core radius of about 14.4pc, eccentricity 0.8 and «=10. 
The sum of these two components is shown in green. The x range 
spanned in the plot corresponds to thrice the value of Xe for this 
model galaxy. 

the raw mass-to-light distribution, while the Figure 9 refers 
to the situation in the smoothed configuration. 

Such a fall-off is also expected for the Sersic galaxies, 
though the exact nature of the trend with an intrinsic scale 
length (such as the core radius in case of the cored galaxies) 
would expectedly differ. But of course, that our recovered 
profiles exhibit this exponential fall-off, indicates the success 
of our modelling. 

7.4 Limitations & Outreach 

The more pertinent question is, how general is this trend 
between the local M/L at the centre and at 3Xe. In other 
words, how convincing is the idea that in general there is a 
quantifiable connection between the luminous matter den- 
sity in elliptical galaxies and the dark matter content at 
approximately 3 times the effective radius? 

The answer obviously depends on the genericness of the 
models that we have consulted in our work. As far as the 
amplitude of the model luminosity densities are concerned, 
it has been constrained to a value that is typical of the order 
of magnitude of the central density of real ellipticals. Also, 
a range of core radii has been scanned successfully in our 
cored models. Though a range of Sersic type model bright- 
ness profiles have not been scanned, the trends noted with 
the single model Sersic galaxy that was considered here, ap- 
pear similar in nature to those found for the cored models. 
Thus, it seems possible that the inner and outer M/L in the 
raw two-stepped M/L profile will be connected even in the 
Sersic case, with the exact nature of this connection given 
by parameters intrinsic to the brightness distribution. 

All this corroborates the generality of our models. At 
the same time it needs to be appreciated that the success 



of the procedure depends on the conformity of a real ellip- 
tical galaxy to the toy models that were used to set up the 
rules that connect the inner and outer mass-to-light ratios. 
In particular, the surface brightness distribution of the real 
galaxy and the compatibility of its dark halo to the NFW 
form need to be checked. In fact for each unique family of 
brightness profiles, a different rule connecting Ti„ and Tout 
is expected to exist. In this paper, we presented the same 
for cored galaxies that were normalised to the same central 
luminosity density and found a similar trend with one toy 
Sersic galaxy. An extensive study of toy Sersic galaxies is 
underway, to form a rule similar to the one prescribed in 
Equation 7. This proposed project will also include an in- 
vestigation of the effect of changing the nature of the dark 
matter distribution in the model galaxies, on this rule. Thus, 
this method can potentially provide a simple automated way 
of deciphering mass profiles of a large sample of elliptical 
galaxies. 

Thus one important conclusion of this paper is that a 
connection exists between the inner luminous matter and 
the dark matter density at about 3i?e in cored ellipticals and 
perhaps also in Sersic galaxies. This notion, though encour- 
aged by intuition, needs to be cultivated better in order to 
ascertain the physical origin of the accompanying parametri- 
sation. 

Similar relationships have been explored in the context 
of disk galaxies; some of more recent works include contribu- 
tions from Sancisi (2004), Verheijen (2001) and Noordermer 
(2006). In these studies and the many more that preceeded 
them, a relationship was noted between the characteristics 
of the rotation curves and the optical luminosities of disk 
galaxies. Investigations along similar lines are of course not 
possible with ellipticals but studies of samples of ellipticals 
in various catalogues have been undertaken to understand 
mass distribution in early type systems (Mamon & Lokas 

2005, Padmanabhan et al. 2004, Lintott, Ferreras & Lahav 

2006, among others). The mass modelling technique that is 
advanced in this paper is one way to pursue the same ques- 
tion, over larger radial ranges than are typically covered in 
these studies (about 3 effective radii as compared to 1). 

Given the simplicity of this method and the twosome 
observational input required, (which are readily provided 
by most catalogues, including the SDSS), a future project 
(Michielsen & Chakrabarty 2006) is planned in which the 
construction of the mass build up in the nearby universe 
will be attempted. 
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